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Abstract 

In this paper we present a unified picture concerning Lie- Trotter 
method for solving a large class of semilinear problems: nonlinear 
Schrodinger, Schroginger-Poisson, Gross-Pitaevskii, etc. This picture 
includes more general schemes such as Strang and Ruth-Yoshida. The 
convergence result is presented in suitable Hilbert spaces related with 
the time regularity of the solution and is based on Lipschitz estimates 
for the nonlinearity. In addition, with extra requirements both on the 
regularity of the initial datum and on the nonlinearity we show the 
linear convergence of the method. 
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1 Introduction 



Let us consider the semilinear evolution equation 



u t + iAu + iB (u) = 0, 
u(0) =w e Hi, 



(l.l) 



where A is a self-adjoint operator in the Hilbert space Hi and B : H x — >■ Hi 
is a locally Lipschitz map. Since a large number of problems fall under this 
situation, at least we can mention the nonlinear Schrodinger, Schroginger- 
Poisson, Gross-Pitaevskii (see [I] for more details), and a large amount of 
articles are devoted to the numerical study of time-splitting methods, most 
of them concerning Lie- Trotter and Strang schemes, see [H El El [71 [8] , we 
shall present in this article a unified picture of time-splitting methods. This 
means that we shall show general results concerning both the order of con- 
vergence, and the regularity required for initial data. Despite the fact that 
we are mainly interested in time discretization, note that the standard result 
for Lie-Trotter schemes developed in the literature expresses that the con- 
vergence is globally linear in the time step, we also take under consideration 
discretization in space (see subsection 13 .4p . In addition, we also show that 
under the (weaker) assumptions made above on the operators the method is 
well defined and converges in the smaller space Hi. 

To see this we first show how to solve the problem (11.11) by means of a 
generic time-splitting scheme. Note that any solution of (11. ip verifies the 
fixed point integral equation 



where <3> A denotes the strongly continuous one-parameter unitary group gen- 
erated by — iA, this means that: v (t) = § A (t) v is the solution of the linear 
problem 



The following well-posedness result of (II. 2p is well-known, for proof and 
details, see [I]. 

Proposition 1.1. Let B be a locally Lipschitz map defined on the Hilbert 
space Hi with B (0) = 0. Then for any u G Hi there exists T* = T*(n ) > 
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and a unique solution u G C ([0, T*(uq)), H x ) of equation (II. 2p . Moreover, 
the map T* : Hi — > C ([0,T*(u )), Hi) zs lower semicontinuous, and for any 
T < T*(uq) the map Hi i— > C ([0,T], Hi) given by wo H- u is continuous, i.e.: 
given e > 0, there exists 5 > such that if \\uq — uq\\ < 5 then T < T* (uq) 
and \\u(t) —u(t)\\ < e for t G [0,T], where u is the solution of (11.21) with 
u (0) = Uq. Finally, is also valid the blow-up alternative: 

1. T*(uq) = oo (u is globally defined). 

2. T*(uq) < oo and lim = oo. 



Since B is a locally Lipschitz map, there exists a flow $ , defined locally 
in time, generated by the problem 



Let $ be the flow of the equation — i(A + B) defined by $ (t) (u ) = u(t), 
where u is the solution of (]1.2p . The idea of time-splitting methods is to 
approximate $, the exact flow, by combining the exact flows Q A and $ B , 
in the following sense: for any (small) time step h > 0, the discrete flow is 
defined by 



where the splitting scheme given by a\, . . . , a m , bi, . . . , b m verifies a± + ■ ■ ■ + 
a>m — bi + - • - + b m = 1. Let us mention that for m = 1 (therefore ai = b\ = 1) 
we get the Lie- Trotter scheme; and for m = 2 and ai — a,2 = 1/2, b\ — 1, b% — 
we get the Strang scheme. Other Yoshida schemes (see details in [TJ]) are 
represented similarly. 

For fixed «o £ Hi and T < T*(uo), the convergence result expresses that 
{u , $/j(m ), . . . , $^(wo)} converges in some sense to the exact solution at time 
t = kh, i.e. {u , $(/i)(wo), • • • , ^( n h)(u )}, when the time step h = T/n goes 
to 0. We note that the splitting scheme given by ai, . . . , a m and b±, . . . ,b m 
is performed n times before reaching the value t = T. Clearly, the scaling 
t — >■ Tt allows us to restrict our attention to the normalized case T = 1, and 
this will be the case in the sequel. We therefore set a, (3 as the 1-periodic 
functions defined by: 



w t + iB (w) = 
w (0) = w . 



(1.4) 



$ ft = $ B (b m h) o $ A (a m h) 



o ■ ■ ■ o $ B (bxh) o $ A ( ai h) , 
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g (t) = !° ,ttj-l<m{t-[t])<j-l/2 
\2mbj , if j - 1/2 < m (t - [t]) < j. 

It is, then, a straightforward computation to verify that for n G N and 
a n (t) = a (nt) , (3 n (t) = (3 (nt) , the continuous flow generated by the (non- 
autonomous) operator — i (a n A + f3 n B), denoted by $„,, verifies $ n (l/n) = 
$/ t . Therefore, the convergence (in time) of the splitting scheme is expressed 
as $ n (t) converges to $ (t) as the time step h = 1/n goes to 0. In what 
follows we shall refer to an abstract time- splitting method when we are given 
a pair of T-periodic functions a, (3. 

Finally, we also take into consideration the convergence in space. It is a 
common practice to solve the problem (11. 3p by means of spectral methods, 
which consists of solving the problem on a finite dimensional invariant sub- 
space (generated by eigenf unctions of the linear operator A). Since invariant 
subspaces of A are not necessarily $ s -invariant, the approximated solution 
is projected before the application of $" 4 ; this gives the (finite dimensional) 
discrete flow: 

$ A = $ B (b s h) o $ A (a s h) o P o • • • o P o $ B (fci/i) o $ A ( ai h) o P, 

where P is the orthogonal projection onto the finite dimensional invariant 
subspace. 

In a more general setting, if we take $^ as an approximation of the exact 
flow this gives the discrete flow: 

$ h = $ B (b s h) o $ A (a s h) o ■ ■ ■ o $ B {bth) o $ A {a x h) . (1.5) 



1.1 Notation and Main Results 

Throughout this paper the evolution problem is given by equation (11.11) 

Ut + iAu + iB (u) = 0, 

U (0) = M , 

for no G Hi, where A is a self-adjoint operator in Hi, and B : Hi — > Hi 
is a locally Lipschitz map. The problem under consideration is to find the 
generated flow $(t) in a compact interval [0,T], where the solution exists. 
The abstract time-splitting method to solve the evolution problem (jl.ip for 
t G [0,T], i.e. to get the flow $(t), is thus described as follows: 
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1. Set a,/3 G L\ oc T-periodic bounded functions with total integral 

/ a= I = 1. 
Jo Jo 

2. Fix n G N and the step size h n = T/n (the choice T = 1 shall be used 
in the sequel). 

3. Set the sequences a n (t) = a(nt) and (3 n (t) = (3{nt). 

4. Get the flow <&h n of the non-autonomous equation u t = —i (a n A + f3 n B) u. 
Under this situation we show: 

Theorem 3.1 (Convergence). Let u G H x and T < T*(u ), then there 
exists no G N such that for any n > n 0; the function $ n (t)w is defined for 
t G [0, T], and lim max II w (t) — u n (t)\\ = 0. 

n->oo tg[0,T] 

In order to get the order of convergence for abstract methods some extra 
regularity both on the time derivative and on the nonlinearity is needed. The 
basic assumption is as follows: let H be a Hilbert space such that Hi C H , 
with continuous embedding, we asume 

1. The solution u of ([I2D verifies u G W 1 ' 00 ([0,T], H ). 

2. There exists a bounded map B' : Hi y B (Ho) such that, for e > and 
u G Hi, the estimate 

\\B (u + w)-B (u) - B' (u) w\\ Ho < e \\w\\ Ha 

holds for some 5 > and for any w G Hy with ||iw|| H < 5. 

Theorem 3.9 (Local error). Let uq G Hi and T < T*(uq), then there exists 
a constant C > and uq G N such that for n > no, the following estimate 
holds for the time step h n = T/n 

||$ (hn) u - $ n (h n ) u \\ Ha < Chi. 

Theorem 3.10 (Global error). Let uq G Hi and T < T*(uo), then there 
exists a constant C > and n G N such that, for n > n : 

max ||$ (kh n ) u - $„ (kh n ) u || H < Ch n . 

0<k<n u 
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2 Auxiliary Results 



This section is devoted to present some basic results that we use to prove 
the convergence theorems. We start with the following notion. We say that 
a sequence {a n } n£ ^ of functions in L\ oc (R) converges weakly to a £ L\ oc (R), 
denoted by a n — ^ a, if for any compact interval 7 C R and 9 £ C (I), the 
following estimate holds 



n— >oo 



lim ^ a n {t) 6{t)dt = J a (t) 9 (t) dt. 



Lemma 2.1. Let a n , a, a £ L\ oc (R), n £ N, sitc/i t/iat a n — ^ a and |a n | < a. 



Then for any 9 £ C ([0, T]) t/ie sequence O n (£) — I a n (£') # (i') dt' converges 



uniformly to 9 (£) = / a (£') 6 1 (£') <it' ; on [0, T]. 



o 



Proof. Suppose n does not converge to uniformly, then there exists e > 
and a subsequence Q n . such that max 10 (t) — n , (t)\ > e. Using the 

0<t<T 

estimate 

|0« fc (*)| < max |0(*)| \\a\\ Ll{[0yT]) , 

we have that the sequence {©n fe } n>1 is uniformly bounded in C([0, T]). 
A similar argument allows us to conclude that the sequence {0n fe } n>1 is 
equicontinuous. By Arzela-Ascoli theorem, we obtain that (a subsequence 
of) Q nk converges uniformly to 0* 7^ on [0, T]. But 0„ fc converges point- 
wise to 0, which is a contradiction. This finishes the proof. □ 

For any real valued function a £ L\ oc (R), we define the propagator op- 

erator § A ' a (t 1 ,t ) = $ A (r (ti, t )), where r (ti, t ) = / a (t) dt. It is clear 



J to 

that the propagator (ti,t ) verifies: 

1. $> A > a (t ,t ) = I. 

2. ^ a (t 2 , t ) = $^ (t 2 , t x ) to). 

3. UueD (A), then d t $ A ' a (t, t ) u = -ia (t) A$ A > a (t, t ) u. 

Observe that if «o £ D(A), then u(t) = <& A,a (t, 0) uq is the solution of 
the linear evolution Cauchy problem iu t = a(t) Au with initial condition 
u (0) = u . 
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Proposition 2.2. Let {a n } n£ -^ be a sequence of real valued functions in 
L} oc (R) such that a n — 1, £/ien (t, t') = $ A ' a ™ (t,f) converges strongly 
to & (t — t 1 ). Moreover, if \a n \ < a E L] 0C (K), then the convergence is 
uniform for t',t on bounded intervals. 

Proof. Let I C R be a compact interval and r n : / x / — > R defined by a n . 
Since a n -± 1, we have r n (£, t') ->■ t—t', thus lim $ A ' n (t, t') u = & A (t - t') u. 

n— >oo 

If \ot n \ < ct, from Lemma [2.11 it follows that the sequence T n (t,t') converges 
to t — t' uniformly on I x I. For any u G D (A), the estimate 

||$ A ' n (t, t') u - $ A (t - t')u\\ < \r n (t, t') -{t- t')\ \\Au\\ , 

is verified. Since D(A) is dense in Hi, using an e/3 argument we finish the 
proof. □ 

Lemma 2.3. Let v G C ([0, T], Hi) and e > 0. Then there exist 8j G 
C ([0,T]) and Zj G H 1; < j < m, such that the function 

z(t) = J2 (2-1) 

0<j<m 

satisfies max \\v (t) — z (t)\\ < e. 

Proof. Let 5 > be such that \\v (t) - v (t')\\ < e/2 if \t - 1'\ < 5, and let 
t_i <to = 0<ti<---<t m = T< t m+ i be a partition with tj — tj_i < 5. 
Let also 8j G C (I) be such that < 6j < 1, Yl @j = 1 an d supp C 

0<j<m 

(tj-i, tj + i). Taking Zj = v (tj) we have for t G [t,_i, t,] 



(t) - z (t)\\ = IK^-i (f) + % (t)) w (f) - 0,-i (t) ^-i - B 3 it) zA\ 

<lk(*)-«(*i-OII + l|u(*)-«(*i)ll- 

Since |t — ^-_i| , |i — tj_i| < 5, the proof is finished. □ 

Corollary 2.4. Let j3 n be a sequence of real valued functions in L\ oc (R) such 
that (3 n ^0 with |/3 n | < (3 G Lj oc (R), and let v E C ([0,T], Hi). Dearie K(t) 
as follows 

V n (t)= [ /3 n (t')v(t')dt' (2.2) 



o 



ThenV n G C([0,T],Hi) and lim max ||K (t) II = 0. 

n->oo te[0,T] 
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Proof. Let e > and let z(t) be the function given by Lemma 1231 We define 

Z n (t)= [ P n {t')z{t')dt' = Q i,n(t) 



Zj, 



0<j<m 



where ir , (t) — / B n (t')9j(t')dt'. From Lemma [27TT lim max \\Z n (t)\\ = 

JO n^oote[0,T] 

0. On the other hand, from Lemma 12.31 we have max \\V n (t) — Z n (t)\\ < 



^ I r 1 1 l 1 ([o r]) wmc h proves the result. □ 

Corollary 2.5. Let v £ C (I, Hi) and {a n } nS N a, sequence of real valued 

functions in L] 0C (W) such that a n — ^ 1 and |a n | < a € L] 0C (M). Then 
$ A,n (t, t') v {t') converges uniformly to <3> A (t — t') v (f) on I x I . 



Proof. Let z(t) be as in Lemma [2.31 then 

($ A > n (t, t') - $ A (t - t')) v (*') = $ A ' n (t, *>) (v (t') - z (If)) 

+ <S> A (t-t>) (v(t')-z(t')) 

+ ($ A ' n (t : t') -$ A (t- t'))z{t') . 

Since $ A ' n (t, t') , $ A (t — t') are unitary operators, the first and the second 
term on the right-hand side are bounded by e. From definition of z, it is easy 
to see that 

<5> A ' n (t,f) -<t> A (t- t')) z (t')|| < max max |0 - (t')\ x 

l<j<m t'&I 



\\{<5> An (t,t') -$ A (t-t')) 

l<j'<m 



Using Proposition I2.2[ we obtain the result. □ 
Let B be a bounded, 1-periodic, function. For n £ N we define f3 n (t) = 

(3 (nt), we note that f3 n — = / /3 (t) dt. Then, under additional hypothe- 

Jo 

ses on v, we obtain an estimate for the order of convergence in Corollary 12.21 

Lemma 2.6. Let v £ W 1 ' 00 ([0, h n ], H), V n {t) be given by Q and w n = 
V n (h n ). If(f3) = orv(0) = 0, thenw n satisfies \w n \ < - \\fi\\ LOO ll^ll z,°°([o,^],H!) K- 
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Proof. Using v (t) = v (0) + / v t (t') dt', we obtain 

Jo 

ph n pt 

w n =(/3)v(0)h n + / / f3 n (t)v t (t')dt'dt 
Jo Jo 

h n pt 

/ f3 n (t)v t (t')dt'dt, 
'o Jo 

ph„ pt 

then ||T4 (h n ) \\ < / / \/3 n (t)\ \\v t (t')\\ dt'dt and an easy estimation implies 
Jo Jo 

the result. □ 



3 Main Results 

3.1 Convergence in Hi 

Let {a n } n€ M, {/9n}nGN be two sequences of real valued functions in L\ oc (M) 
such that a n , (3 n — ^ 1, \a n \ < a and \(3 n \ < (3, with a, (3 G Lj" oc (1R). For n G N 
we consider the approximated evolution problem, 

iw t + (a n A + (3 n B) w = 
w(0) = u 

related with the abstract splitting scheme defined by these sequences, and 
we denote by $ n the related flow. (The exact flow will be denoted by $.) 
Let Mo G Hi be given and let u n = $ n Mo be the solution of the problem (13.11) . 
we recall below the integral expression for u n 

u n (t) = $ A ' n (t, 0) uo - i I n (f) <$> A ' n (t, t') B (u n (*')) dt' . (3.2) 

Jo 

We are now in position to give the first result concerning the uniform 
convergence of $ n (t)tio to $(t)iio for t G [0,T] and for any uq G Hi. 

Theorem 3.1 (Convergence). Let Uq G Hi and T < T*(uq), then there 
exists no € N such that for any n > uq, the function § n (t)ito is defined for 
t£ [0,T], and lim max ||« (t) — u n =0. 

n-voo f€[0,T] 

Proof. For t < min {T, T*(uq)}, we write 

u (t) - u n (t) =h tn (t) -i ($ A (t) J 2 , n (t) + J 3 , n (t) + J 4 ,„ (*)) , (3.3) 
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where 

I 1>n (t) = ($ A (t)-$^ n (t,0))« , 

/ 2 ,n (t) = T (1 " & (0) (-*') 5 (« (*')) dt', 
JO 

h,n (t) = f n (f) ($ A (t - f) - (*, 0) B (u (f )) < 

/4,» (*) = / /?n (*') $ A ' n (*, (B (u (0) - 5 K (0)) dt'. 



We shall prove that ij >n (£) — y as n — > oo uniformly on [0, T\. 

From Proposition 12.21 we get lim max ||/ ln (t)|| =0. From Corollary [2^1 

n-s>oo te[0,T] 

we deduce lim max ||/2n(£)ll — 0. 

n-Kx>ie[0,T] 

For j = 3 we have the estimate 

H J 3,n Wll < PL ([ 0,T]) ^ II (* - " (*> 0) B (U (0)|| • 

Using Corollary 12.51 we obtain lim ||/ 3n (t)|| = 0. 

Let R = max II w and let L be some Lipschitz constant of B on the 

te[o,T] 

ball of radius 2R centered at the origin. Then there exists n G N such that 
for n > n is valid the estimate 



max 

te[o,T] . 

3=1 



X] Il4n(*)|| < £ex P (~ L PLi 



([0,T]) 



Thus, we have 

\\u (t) - u n (t)\\ <eexp y—L 



Li([0,T]) 



+ L / 0(0 ||u(0 -u« (OIK 
Jo 

from Gronwall inequality we obtain \\u (t) — u n (t)|| < e, and then T < 
T*(u ). This finishes the proof. □ 



3.2 Error estimate 

In this section we obtain local and global in time error estimates for general 
time-splitting methods. These results are optimal for Lie- Trotter schemes, 
whose local convergence in the whole space is quadratic in the time step. 
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Let a, /3 be 1-periodic, bounded functions, with (a) = (f3) = 1, and set 
a n (t) = a (nt), (3 n (t) = (3 (nt), with h — 1/n \, 0. We recall that, under this 
situation a n , (3 n — ^ 1. In order to get these error estimates we impose some 
regularity both on the time derivative of the solution and on the nonlinearity 
B, which is accomplished as follows. We consider a Hilbert space Ho such that 
Hi is continuously embedded in H , and there exists a self-adjoint extension 
of the operator A : D — > H with H x C D. We can see that for uq G H 1; 
the solution u of (JH2D or fl32D verifies u G W 1 ' 00 ([0,T], H ). We also assume 
that there exists a map B' : Hi i— > B(Hq) such that for R, e > 0, it can be 
chosen C, 5 > verifying 

\\B'(u)\\ B{Ho) <C, (3.4a) 
\\B{u + w)-B (u) - B' (u)w\\ Ho < e |MIh„ . ( 3 - 4b ) 

for u,w G Hi, \\u\\ H < R and ||w|| Hl < o~- 

From conditions (13.41) it is clear that for R > 0, there exists L > such 
that \\B (u) — B (v) || Ho < L \\u — v || Hq for any u, v G Hi with ||w|| Hl , \\v || Hi < 
R. Let w , u G Hi, T < min {T* (u ) , T* (u )}, e > and R > such that 

||$ (t) wo|li,«»([o ) r| j H 1 ) ' ll $ (*) «oIIl«»([o,ti ) Hi) - R 
since $ A (t — t') {Q A ' n (t, t')) is an unitary operator of H , we deduce that 

||$ (t) U - $ (t) M || Ho < \\ U ~ M || Ho + L [ II $ (0 MO - $ (0 «0 || Ho 

Jo 

Therefore, we have the estimate 

||$ (t) u - $ (t) wo||h < e L < ||«o - ^oIIho • (3-5) 

We now define for a fixed T > the space X T = C ([0, T], H^rW 1 ' 00 ([0, T], H ). 
Since £> is a locally Lipschitz map and conditions (13.41) we can see that 
u i — y B o u is Sl well-defined bounded map in X T and (B o u) t = B' (u) u t . 

The following lemma deals with local nonlinearities. 

Lemma 3.2 (Local nonlinearities). Let f : C — > C be a smooth map in the 
real sense, (i.e.: iff = then the map (£,77) H- (/« (f + mj) , /« (£ + z't?)) 

is smooth on R 2 ). Let also H x = H s (R d ), with s > d/2, and H = L 2 (R d ). 
Then B : Hi 1— > Hi given 61/ -B(w) = /(«) is « well-defined map, in addition 
B' : Hi 1— > B (Ho) given by B'[u)(v) = f'(u)v is well-defined and verifies 

rn- 



Proof. From Schauder lemma (see Theorem 6.1 in [H]), for s > d/2, it follows 
that B : H s (R d ) \-> H s (R d ) is a well-defined, locally Lipschitz map. Taking 
norm in the identity 

/' («) = ( f { : ] (u) + /« («) + < Of («) «;« + f? («) ™« 



<C Hix 



, with C(#) = max |f 



we obtain («) w|| £ 2( R d) S o (J|u|| £ oc( Rd ) J ||w|| £a 

Using |/ (u + w) — f (u) — /' (u) .w\ < e \w\ if \u\ < R and \w\ < 5, we get 
the required inequality. This finishes the proof. □ 



In order to add Hartree-type nonlinearities we first collect some useful 
estimates. 

Lemma 3.3. Let W x G L°° (R d ) ,W 2 G LP (R d ) , with p > 2, p > d/A. 
Let also u G H s (R d ) , with s > d/2, and v G L 2 (R d ) . Then the following 
estimates do hold, with C depending only on s: 



(i) ll^i* Re (Vt; 

(%%) ||H/ 2 *Re(w*w 

(Hi) || Wi * \u\ 2 \\ Loo 

(iv) \\W 2 * |u' 2 " 



Proof. Estimates (i 



| £ oo( R d) < ||W^l|| LO o( R d) IMI £ 2( R 

\ L co( R d\ < C\\W 2 \\ LP ( ud ) \\v\\ L2( 



u 



H"(R d ) 



u 



Rd )<C\\W 2 \\ 



u 



Li' l 



U 



20 
\H' 



I 2(1-0) 



and (iii) follows immediately from Young and Holder 
inequalities, while estimates (ii) and (iv) also uses Gagliardo-Nirenberg in- 
equality. □ 

Lemma 3.4 (Hartree-type nonlinearities). Let W G L°° (R d ) +L P (R d ) , with 
p>2,p> d/A, let Hi = H" (R d ), with s > d/2, and H = L 2 (R d ) . Then 
B : Hi y Hi, with B(u) = (W * \u\ 2 ) u is a well-defined map, in addition the 
map B' : H x ^ B (H ) given by B\u){v) = (W * \u\ 2 ) v + 2 (W * Re (u*v)) u 
is well-defined and verifies estimate fjOj. 



Proof. Since 

B(u + v) - B{u) =(W* \u\ 2 ) v + 2 (W * Re {u*v)) u 

+ 2 (W * Re (u*v)) v + (W * \v\ 2 ) (u + v), 

the linear term is given by B'(u)(v) = (W* \u\ 2 ) v + 2(W *Re(u*v))u. The 
estimate ( 13 .4p follows directly from Lemma 13.31 □ 
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Theorem 3.5 (Local error). Let u G and T < T*(u ), then there exists 
a constant C > and n G N such that for n > uq, the following estimate 
holds for the time step h n = T/n 

||$ (h n ) u - $„ (h n ) u \\ Ho < Ch 2 n . 

Proof. Replacing t = h n in Eq. (13. 3p and using that $ A (/i n ) are unitary 
operators, we see that it is sufficient to show the estimates \\Ij in (h n )\\ H < 
Ch 2 n , where Ij >n are defined as in Theorem 13. 11 Since (a) = 1, we have 

h,n (K) = ® A (K) - (h n , 0) = $ A (h n ) - ® A (h n {a)) = 0. 

From Theorem 13.11 there exists Uq G N such that for n > n it holds 

T* > T and max < max \\u(t)\\ + 1 = R. Setting r/ 2 > (i) = 

te[o,T] te[o,T] 

$ A (-i) 5 (u (t)), it is clear that v {2) G X T and 

vf ] (t) = $ A (-t) (iAB(u(t)) + (B(u(t))) t ) : 



from where it follows the estimate 
Using that 



,(2) 



L°°([0,fe„],H Q ) 



< C(i?). 



/ ■ -ft 

I 2 ,nW= / (l-/3„(s))^( S )d S 
./o 

and since (1 — /3) — 0, from Lemma [2.61 we deduce 

\\l2,n(h n )\\ Ho < C (R) (1 + \\P\\ LOO ) hi 

We set (t) = ($ A (h n - t) - $ A - n (h n , t)) B (u (t)). It is clear that G 
X T , w (3) (0) = 0, and 



,(3) 



(t) =i ($ A (K -t)-a n (t) § A ' n (h n , t)) AB (u (t)) 

+ ($ A (h n - 1) - $ A <" (h n , t)) (B (u (t))) t . 



,(3) 



L°°([0,h n ],H ) 



<C(i?)(l 



Taking norms, we deduce the estimate 
Using Lemma [2.61 again, we obtain 

\\I 3 ,n(h n )\\ HQ <C(R)(l + \\a\\ L „)\\f3\\ L „h 2 n . 

We finally set v& (t) = $ A ' n (h n , t) (B (u (t)) - B (u n (t))). Since 

(t) = ia n (t) <$> A < n (h n , t) A (B (u (t)) - B (u n (t))) 
+ $ A > n (h n ,t) (B (u (t))-(B (u n (t)))) t 

and u, u n are bounded in Xt, using a similar argument as in previous cases 
we deduce the estimate for (h n ). Theorem is thus proven. □ 
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Under hypotheses of Theorem 13.51 we formulate the result concerning 
global error estimate. 

Theorem 3.6 (Global error). Let uq G Hi and T < T*(uq), then there exists 
a constant C > and Uq G N such that, for n > n : 

max ||<I> (kh n ) u - $„ (kh n ) u \L < Ch n . 

0<k<n "° 

Proof. Setting e k = ||$ (kh n ) u — $ n (kh n ) Mo|| H , it follows that 
e k +i < ||$ (h n ) $ (kh n ) u - $ (h n ) $ n (kh n) ^0 1| Hn ' 

||$ (hn) ®n (kh n ) U - $ n (h n ) ($ n (kh n) u 0) || H () • 

Using estimate (13. 5p and Theorem I3.5[ we deduce e k +i < e Lhn ek + Ch^, from 
where, by means of an inductive argument, we conclude the estimate, valid 
for < k < n, 



fc-i 

e k <Ch 2 n Y,e L3hn 

3=0 

This finishes the proof. 



Cfl n ( p Lkh n _- [ \ < C ( eLT - l ) , 



□ 



Corollary 3.7. Let \-\ e = [H , Hjjg be the interpolation Hilbert space, 6 G 
(0, 1) and Uq G H . IfT<T* and s > 0, then there exists n Q G N such that 

max \\$(kh n )u - § n (kh n )u \\ H < e h l ~ e , 
holds for n > n Q . 

Remark 3.8. Let u ,u G H , and let T < min{T* (u ) ,T* (u )}. Using the 
notation and the result of Theorem 13.61 and the estimate (13. 5p we deduce 

||$ (kh n )u - $ n (kh n )u \\ Ho < ||$ (kh n )u - $ (kh n )u \\ Ho + 

||$ (khn) U - $„, (khn) M ||h 



<e LT \\u - uq\ 



h + Ch n . 



3.3 Approximation methods 

Assume we can define an approximation $^ for the flow $ A such that for 



any u G Hi, 
step h, 



Q A (t)u < C ||w|| Hl and for any uq G Hi and a small time 



$ A (h)uo-$ A (h)u <C/i 2 ||t/o|| Hl - 

Ho 



(3.6) 
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Let $ ft be the flow given by (11.51) . From the identity $ A (t) = $ A (t) + 
(t) — $ A (i)Y we get the following decomposition for the discrete flow: 

$ ft = $ h + N h , where 



7 e{o,i} s i=l 
7^0 



Proposition 3.9 (Approximation method). Let $ be an approximation of 
the flow $ A satisfying (13.61) . Lei uq G Hi, T < T*(w ) ; #ien t/jere exists a 
constant C > and no € N suc/i t/iat 7 /or n > no '■ 



max 

0<fc<n 



$ (kh n ) u - $*u 



H 



<Ch n . 



Proof. Using that B : Hi — > Hi is Lipchitz with constant L, then for all 
u G Hi and for all s 



\<S> B (b s h)u\\ Hi < e 



L{b s h) 



H, ' 



which combined with inequality (13.61) yields ||-WfcUo||H < Ce Lh h 2 \\uq\\ h . Us- 
ing that 



$ (h) u - <5> h u < ||$ (h) u - $hU \\ H 

Ho 



$^0 - $hUQ 



H„ 



and theorem (13. 5p . we obtain that there exist n such that for n > n 

$ (h n ) u - $ hn u 



H„ 



<Chl 



and therefore we deduce the desired inequality. 



□ 



3.4 Spectral methods 

We then turn to the discretization in space variables. Let R > be fixed, let 
E be the projection valued spectral measure of A : H x C D (A) — > H , and let 
P = E([—R, R]) be the orthogonal projection onto the yl-invariant subspace 
H = P(Ho). According to previous subsection, we define $ A = & A o P and 
$ a ^ = $ a ^ ( P + j _ P ) = $a ^ + $ a q (I -P). We get the following 

decomposition for the discrete flow: $^ = $^ + A^, where h > is a small 
time step and 

s 

jVfc= Yl Y[® B (fyh) o $ A ( aj h) p 1 -^ (i - py s . 

7 e{o,i} s j=l 
7^0 
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Theorem 3.10 (Spectral approximation). Let u G H 1; T < T*(uq), and 
n G N be given. Then, for R > h~ 2 = (n/T) 2 is valid the estimate: 



max 

0<fc<n 



< Ch r 

Ho 



$ (kh n ) u - $> 
Proof. For any u G Hi we have 

||u - Pu\\l o = [ d(u\E (A) u) Ho < R- 2 [ \ 2 d (u\E (A) u) Hq 

J\X\>R J\X\>R 

and then \\u — Pu\\ H < R^ 1 \\u\\ H . Being $ A a unitary operator, we get that 
— P)u|| < R^ 1 ||u|| Hl . Taking R > h~ 2 we get the desired inequality 
from proposition (13. 9p . 

□ 

When (A ± z) -1 are compact operators, there exists a basis {^} >0 C 
D (A) of H and a sequence {Ay} >0 C R with |Ay| f oo such that A<£>j = A^.,-. 
The operator <3> A (£) P could be written as 

$ A (t)Pu= ^(vMhoW' 

which represents the approximate solution of (11. 3p in terms of the eigenfunc- 
tions (which in most cases are explicitly given). 



4 Examples 

4.1 Nonlinear Schrodinger equation 

We consider 

iu t + Au + f{\u\ 2 )u + (W(x) * \u\ 2 ) u = 0, 

U (0) = M , 

where / : C — > C is smooth as a real function, and W(x) is an even func- 
tion such that Wi G L°° (R d ) ,W 2 G L p with p > 2, p > d/4. Taking 
Hi = H s (R d ) and H = L 2 (R d ), with s > d/2, s > 2, we can see that 
A = —A is a self-adjoint operator, and B (u) = — f(\u\ 2 )u — (W(x) * |w| 2 ) u 
is a locally Lipschitz map (see, Lemmas 13.21 and 13.31) . Following these lem- 
mas we can also deduce that, for any uq G H 1; and T < T*(u ), the solution 
verifies u G W 1,oc ([0,T], Ho) ; in addition, the nonlinearity B satisfies (13.41) . 
We thus obtain Theorem 4.1 of [3] for Lie- Trotter splitting schemes. Us- 
ing H 2e (R d ) L°° (R d ) for 6 > d/4 and Corollary [321 we can see that 
\\u (kh) - u ri {kh)\\ Loo r Rd \ = o (h 1 ' 9 ). 



16 



Remark 4.1. Since, for d = 3, the Newtonian potential W(x) = \x\ _1 verifies 
the hypotheses of Lemma 13.3} the convergence results are also valid for the 
3-D Schrodinger-Poisson equation: 

iu t + Am + Vu = 
AV = -\u\ 2 

Remark 4.2. In lower dimensions, d — 1,2, the kernel W is not bounded and 
therefore Lemma 13.31 does not apply. Actually, the existence of dynamics 
requires some extra work, see O [10], mainly connected with a suitable de- 
composition of the nonlinearity. However, the conclusions of Theorem 13.11 
13.61 remain valid but their proofs are more involved. 

4.2 Gross-Pitaevskii equation with a trapping poten- 
tial 

We consider the d- dimensional initial value problem 

iu t + An — Vtu — |n| 2 n = 0, 
n (0) = w , 

where Q is a positive definite quadratic form. Without loss of generality we 
can assume Q (x) = u\x\ + - ■ - + u^x d . This equation is used to describe Bose- 
Einstein condensates. The operator A = — A+Q has a basis of eigenf unctions 
(explicitly) given by 

d 

d 

for k = (ki, . . . ,kd) G Nq with eigenvalues Ak = d + 2 ^ ^jOJ 2 , where <fk 

3=1 

is the k-th Hermite function. In [7J the convergence of a split-step method 
using Hermite expansion is studied, the Hilbert spaces H s (lR d ) = D (/W 2 ) 
are defined as the functions u in L 2 (lR d ) such that ||w||^/ Rd \ is finite, where 

Since A > -A, we see H 2 (M 2 ) ^ H 2 (R d ), in particular H 2 (M. d ) 
L°° (M d ) if d < 3. In these cases, Lemma 2 in [7] implies D (A) = H 2 (M 3 ) 
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is an algebra and then B (u) = \u\ 2 u is a locally Lipschitz map. Using 
similar arguments as in the proof of Lemma 13.2} we get (13.4)) for the cubic 
nonlinearity. Therefore, taking Hi = H 2 (M 3 ) and Ho = L 2 (M 3 ), we obtain 
the convergence result given by Theorem 13.61 and like in the example above 
||u (kh) - u n {kh)\\ Loc ^ = o (h e ), for < 1 - d/A. 

Lemma 4.3. For any u G D (A) the following estimate do hold: 

c" 1 (Au\Au) L2(Rd ^ < \\-Au\\ 2 L2 ^ + ||ftu|£ a ( Rd ) <c{Au\Au) L2 ^ 

with c = max {2, 1 + 2d~ 2 Y? j=1 

Proof. Since S (M d ) is dense in D (A), we just have to prove the norm equiv- 
alence for any Schwartz function 

( Au \ Au ) L2(Rd) = \\- Au \\%(Rd) + IMlJ^Rrf) ~ 2 ( Au P U ) L ^(r^ 

Using (Au\Qu) L2 ^ Rd ^ = — (VQ.'Vu\u) — (f2Vii|Vii), we get 

2 (Au\Qu) L2(Rd) < - 2 (VQ.Vu\u) L2(Rd) = (AQu\u) L2(Rd) 



d 

2 



U W L 2( 



Since (Au\Au) L2 f R d\ > d 2 ||u|| i2 ( R d), we have 

||-Au|| 2 2(Rd) + \\nu\\ 2 L2{Rd) < (l + 2cT 2 £ d j=1 ufj (Au\Au) 
From 2 (Au\Q.u) i2 / Rd \ < ||-Aw||^ 2 / R ^ + \\nu\\\ 2 ^ Rd y we obtain 

(Am|Am) l2 ( k ^ < 2 H-AuHja^ + 2 ||nu||^,^ 

and then, the lemma follows. □ 
Corollary 4.4. For d<3, H 2 (R d ) is an algebra with the pointwise product. 
Proof. From the estimate ||fhM;|| i2 ^ Rd \ < ||fiw|| L2 ^ |M| LOO ( Rd ) and the em- 
bedding^ 2 (R d ) L°° (R d ) , we obtain \\ttuv\\ L2 ^ < C ||u||w r< a |M|£ 2 ( Rd y 
Using —A (uv) = —Au v — uAv — 2Vu. Vw, we have 

H- A (HII L 2(Rd) < l|- Au lli2(H<J) IK ; ll L -(Rrf) + \\- Av \\ L 2(R*) \\ U \\l°°(**) 
+ 2 \\^ U \\ L ^R^ ||Vu|| L 4( K d) • 
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Since 



I 2 sr> iu,ii( 4 ~ d )/ 4 II A..n( 4 + d )/ 4 



||V M || i4(Rd) <C\\u\\ ( >< ||-Ati|, ( . 

/ 2 2 \ 2 

< C (JMIl^r*) + \\- Au Wl^) ) ^ C W u Wh*(w<i) » 
we get ||-A(W)|| L 2( Rd ) <C\\u\\ &2 ^ \\v\\ 62 ^ and 

IMIIi^Rd) - C H n lli? 2 (R d ) IMI#2(r<*) ' 

this finishes the proof. □ 



Proposition 4.5. Let f be as in example \3.S\ and d < 3, then the map 
u i— > / (u) is bounded and locally Lipschitz on H 2 . 

Proof. Let R > such that ||it||jrW K <n < R, since \f (u)\ < C\u\ if \u\ < R, 
we have (u)\\ L2 ^ < C \\Qu\\ L2 ^ Rd y Using that Af (u) = f" (u) |Vw| 2 + 
/' (u) Au, we obtain 



\\-Af(u)\\l 2{Rd) + \\nf(u)\\ 2 L2M <C (\\-Au\\ 2 L2{ 

+ \\Qur L2{Rd) + \\Vu\\ 

L 4 (R d ) 

from ( 14.1 p and Lemma 14.31 we have 

(Af (it) |A/ (u)) L2{Rd) <C \\-Af {u)\\ 2 L2(Rd) + {u)\\ 2 L2(Rd) 

<C (\\u\\ 2 L2[Rd) + ||-A U ||^ (Kd) ) < C (Au\Au) 

Let u,v E H 2 (R d ) such that ||«||^2^ , H M llH2( K d) < R, then 



\\f (u) - f (v)\\ &2 ^ < / ||/' ((1 -t)u + tv)\\ m(Rd ^ \\u-v\\ &2 ^ 



dt 



<C\\u-v\\ &2 ^ 



which expresses that / is a locally Lipschitz map. □ 

Using similar arguments as those used in the proof of Lemma I3~2| we can 
see that the nonlinear local term given by B (u) = f(\u\ 2 )u verifies (13.41) and 
then the conclusion of Theorem 13.61 holds. 
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4.3 Nonlinear wave interaction model 



Consider the system of evolution equations modelling wave-wave interaction 
in quadratic nonlinear media (see [2] and references therein). This model 
describes the nonlinear and nonlocal cross-interaction of two waves in 1 + 1 
dimensions. The interaction is described by nonlocal (integral) expressions: 

vgu^ = 0, 
- ug * u W = 0, 

4 1} ,u (2) (0) = 4 2) , 

where v — ±1 and g x = u^>* g(%) — > when x — > — oo. Consider 
the spaces Hi = H 1 (R) x H 1 (R), H = L 2 (R) x L 2 (R) and the operator 
A = id x a z . Define B {u) = vg [u] a y .u, with 




g(u)(x,t) = / U W*(y,t)uU(y,t)dy 
and <J y ,a z the Pauli matrices. Taking 

(g'(u)w)(x,t)= [ X (w®* (y,t)uU (y,t) (y,t)wW (y,t)) dy, 



we can see that B' (u)w = vg' (u) w a y .u + vg (u) a y .w. From Cauchy in- 
equality, we get \\g' (u) w\\ L<x , R -, < \\u\\ L2 ^ HHIz,2(kv From the expression of 

B'(u)w, we conclude \\B' (u) w\\ L2 ^ < C \\u\\ 2 L2 ^ HHIl 2 (r)- Then, (13.41) is 
verified and therefore the conclusions of Theorem 13.51 and Theorem 13.61 are 
valid. 

As an application of these results, we study the behavior of solutions 
with compact support. If supp (wo) C (a,b), since A is a first order lin- 
ear wave equations and it holds supp (£> (w)) C supp(w), it follows that 
supp (<£ A (t) uq) C (a — t, b + t) and supp (<3> B (t) u) C supp (u). Therefore, 
supp {u n {t)) C (a — t, b + t) which implies supp {u (t)) C (a — t, b + t). 

5 Numerical example 

Consider de Schrodinger-Poisson equation in T, i.e. u is a 1 -periodic solution 
of 

iu t + u xx + \u\ 2 u + Vu = 0, 

V xx = V-\u\\ (5.1) 

U (0) = Uq, 
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where T> e C°° (T) is a given real-valued function. We assume that neutrality 
condition is verified: 



/ V (x) dx = \\u \ 
Jt 



2 

L 2 (T) ' 



since \\u (t)||^2( T ) is a conserved quantity, this condition holds for any t. The 

potential V can be calculated by V = —G * g, where g = V — \u\ and G 
is the Green potential defined as the 1 -periodic function such that G (x) = 
x (1 - x) /2 on [0, 1]. We consider H = L 2 (T), Hi = H 2 (T), defining the 
self-adjoint operator A = —d xx and 

B (u) = — \u\ 2 u + (G * g) u, 



we can write (15. 1H in the form (II. ip and from Lemma 13.31 B verifies (13.41) . 
Th 

where 



The linear flow $ A can be written as (t) u) (x) = J2 P & ^p e ^ P t e t2lTpx , 



u p = I u (x) e~ t2jTpx dx. 

T 



Let w be the solution of (11.41) with w (0) = u, using V is a real-valued 
potential, we can see that Re (w*w t ) = 0, which implies \w\ = \u\ and then 

V is constant in t. Therefore <3> s [t)u = e l< ( y+ ' n ' )u, where V is calculated 

2 

using u. Observe that if g = T> — \u\ , then it holds ^ — and the potential 
can be expanded by V (x) = — J2 P ez Gp (^^p) -2 e l2npx . 



5.1 Solving by Discrete Fourier Transform 

We show a numerical method using discrete Fourier coefficients. Let m be 

1 * 

the odd integer m = 21 + 1 and consider (I m u) (x) = U p e' l2 ' Kpx , where U p 

P =-i 

is the discrete Fourier coefficient given by 



1 m— 1 

U p = — f/ e ~ l27rp<?/m 



g=0 



and C/ 9 = u (g/m). Since e " i27rM / m = e - i2 ^P ±m V m , we have f> p = U p±m . We 
also know that 



m— 1 

Tpq/m 



p=0 



It is known that ||u — J m u|| L2 ( T ) < Cm 2 ||w||#2( T ) (see Lemma 2.2 in [12] 
and then we have 



21 



Proposition 5.1. Let $ ? „ (t) = <J> (t) I m , for any u G H (T) it is verified 
\\$ A (t)u- $i (t) u|| x2(T) < CrxT 2 ||«|| H2(T) . 

We can see (t) is an approximation of the flow <3> A that verifies inequality 
( 13.61) in subsection 13.31 for m > n. From definition of $^ (t) and £/ p = U p ± m , 
it holds 

-i47r 2 p 2 t i2npq/m 



p=-i 

m—1 £ 
^ ^ jy- g— i47r 2 (m— p) 2 1 ^i2irpq / m _|_ ^ ^ ^—i4Tr 2 p 2 t^i2Trpq/m 

p=l+l p=0 
m—1 

p=0 



-iXpt i2irpq/m 



where A p = 4m 2 n 2 h (p/m) for < p < m — 1 and h (u) = v 2 — 2 (v — 1/2) 
The solution of (11.41) can be exactly calculated as 

($ B (t) u) (q/m) = e it{v " +N ^U q , 
where N q = \U q \ and the potential V is given by 

m—1 

v q = - Q P K 2el2np9/m > 

p=l 

with q p = Dp — N p . Observe that the neutrality condition reads as go : 
V — N = 0. Therefore, the Lie-Trotter algorithm can be written as: 

- Fix n. 

- Asign h — Tjn. 

- Fix m ~ h^ 1 . 

- Transform D to D using FFT. 

- Compute A -2 . 

- Compute exp (—i\h) . 

- Evaluate U = Uq (q/m) for q = 0, . . . , m — 1 . 
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- For k — 1, . . . , n do 



1 . Transform U to U using FFT (m x log (m) ops) . 

2. Multiple U by exp (— iXh) (m ops). 

3. Obtain £/ (A) anti-transf orming FFT e~ iXh .U (m x log (m) ops) 

4. Compute iV = |[/ (A) | 2 (m ops). 

5. Transform N to N using FFT (m x log (m) ops). 

6. Compute g substracting N from D. 

7. Multiple q by A -2 (m ops). 

8. Obtain V anti-transf orming FFT — \~ 2 .g (m x log (m) ops). 

9. Sum iV and 1/. 

10. Evaluate exp (ih (V + N)) (b x m ops). 

11. Obtain U multiplying exp (ih (V + N)) .U^ A) (m ops). 

12. Asign C/[fc] = U. 

The computational cost is proportional to n x m x log (m). 

To illustrate Theorem (13.61) we present a numerical experiment in one 
space dimension. We use the algorithm described above to discretize the 
Schrodinger-Poisson equation (15 .ip with initial data uq(x) = sin2 +Q (7rx) with 
a > small so that m G H 2 but «o ^ H 2+s for s > a, and X> (x) = 
7 (a) (1 + (1 + 16tt 2 )) cos (4vrx), with 

T(a + 2) 

7(a) ~ 



v^r(a + |)' 

Figure 1 shows the order dependence of the L°° error at time T = 1 on the 
time step-size h. The calculations are performed with a space discretization 
of 2 x 10 5 + 1 and compared to the result with a time step-size h = . 

Figure 2 illustrates the dependence of the L°° error on the space dis- 
cretization parameter n. Here, we use a fixed time step-size h = 10~ 3 and 
compare the results with the result for n = 2 14 + 1. 
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Figure 1: Discretization error for different time step 
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Figure 2: Discretization error for different space discretization 
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